Relativistic Contributions to Deuteron Photodisintegration 
in the Bethe— Salpeter Formalism 



K. Yu. Kazakov* and S. Eh. Shirmovskyt 
Laboratory of Theoretical Nuclear Physics, 
Far Eastern State University, 
Sukhanov Str. 8, Vladivostok, 690600, Russia 
(February 8, 2008) 



Abstract 



In plane wave one-body approximation the reaction of deuteron photodis- 
integration is considered in the framework of the Bethe-Salpeter formalism for 
two-nucleon system. Results are obtained for deuteron vertex function, which 
is the solution of the homogeneous Bethe-Salpeter equation with a multi-rank 
separable interaction kernel, with a given analytical form. A comparison is 
presented with predictions of non-relativistic, quasipotential approaches and 
the equal time approximation. It is shown that important contributions come 
from the boost in the arguments of the initial state vertex function and the 
boost on the relative energy in the one-particle propagator due to recoil. 
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I. INTRODUCTION 



Recently experiments on elastic and quasi-elastic lepton-nucleus scattering have been 
considered as one of the fundamental and reliable resources of probing structure of nuclei 
over a wide range of energy-momentum transfer. A specific place among such reactions is 
taken by deuteron photodisintegration. The reaction 7 + d — > p + n has been thoroughly 
investigated in the region of small and medium energies of the incoming 7-quantum (an 
excellent review of experiments and theoretical frameworks applied to the study of the 
process can be found in Ref. IJ). The theoretical analysis of experimental data obtained 
in these kinematic regions produced an important information on deuteron wave function, 
allowed one to discriminate various contributions to the differential cross section which 
stem from meson-exchange currents, isobar configurations, spin-orbit and further relativistic 
corrections. 

At the present time deuteron photodisintegration is rated among the leading trends at 
experimental facilities around the world. Here it is worth-while to point out a performed 
experiment with the polarized LADON gamma ray beam to investigate the existence of 
narrow dibaryonic resonances in reaction on nuclei at a low excitation energies, an experimen- 
tal program carried out on the linearly polarized photon beam at YEREVAN Synchrotron to 
study the cross section asymmetry of the deuteron photodisintegration process in the energy 
range 0.9-1.7 GeV for proton center-of-mass angle of 90° ||. This is connected the problem 
of the validity of the constituent quark counting rules at energies of a few GeV. Present 
experiments at SLAC with photons of an energy 2.8 GeV [|J and proposed experimental 
programs in RCNP Cyclotron with photon beam at an energy up to 8 GeV [j|] allow one to 
focus on investigation of hadronic systems at quark level. 

Future studies at the TJNAF are approved to extend measurements on the differential 
cross section to a wide range of reaction angles at high energies. The first measurements of 
the cross section from 7 + d^p + nupto 4.0 GeV are in good agreement with previous 
low energy measurements [H] . A comparison of the high precision Mainz data ]7||| (over the 
photon energy range 100-800 MeV) with meson-exchange models incorporating relativistic 
effects P-|TT|] shows that one still lacks of the ultimate conclusion on the role of non-nucleonic 
degrees of freedom in nuclei. It is also understood that further development of deuteron 
photodisintegration theory is needed and consistent treatment of relativistic effects have 
to be applied. This demands a construction of a genuine relativistic formalism for the 
description both the deuteron structure and the reaction mechanism. 

The formulation of a completely relativistic formalism of hadronic bound states and reac- 
tion with them can be developed on the basis field-theoretical Bethe-Salpeter (BS) equation 
for the nucleon-nucleon (NN) scattering [I2|. However approximate methods evolved from 



the BS formalism due to substantial mathematical and computational problems. We bear 
in mind the quasipotential (QP) approach, which reduces the 4-dimensional BS equation to 
a relativistic 3-dimensional equation: the Logunov-Tavkhelidze [DJ, Blankenbeckler-Sugar 
(BbS) JTJJ] and Gross |T3[ and other approximations. 

Actually the relativistic description of the reaction with the deuteron is extensively de- 
veloped in the framework of the BS formalism, which is explicitly Lorentz covariant, provides 
two-body unitarity and includes nucleon and anti-nucleon degrees of freedom in a hadronic 
state in a symmetrical manner. So far applicability of the BS equation to reactions with the 
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deuteron has been bound to: elastic electron-deuteron scattering \12\, elastic pd-backward 



scattering ||16| , inclusive quasi-elastic electron-deuteron scattering |T7|J18[1 and description of 
the static properties of the deuteron (p^EO]. The common feature of these processes is that 
the reaction amplitude in the impulse approximation is proportional to the averaging of cur- 
rent operator between deuteron states. On the contrary, the reaction amplitude for deuteron 
photodisintegration is a non-diagonal matrix element between the incoming deuteron and 
outgoing 2N state. 

An approach towards a covariant description of the reaction 7 + d — > p + n, in which 
the basic degrees of freedom are taken to be hadronic, is developed in the framework of 
the dispersion relation technique for laboratory photon energies E 1 < 400 MeV [J2TJ] (this 
technique is appropriate for the analysis of partial wave amplitudes and takes into account 
final state rescatterings) . Other approach describes deuteron photodisintegration as a simple 
parameterization of covariant deuteron in terms of a hard component and imposing gauge 
invariance on the cross section [^] . This approach looks into the onset of scaling in exclusive 
photodisintegration for high energies in the range 1-4 GeV, as it can be brought about by 
mechanisms which are different from those of pQCD. 

The BS formalism presents a separate view on the problem. Formal application of the BS 
formalism to deuteron electro- and photodisintegration can be found in Ref. |23[| . Here the 
rigorous derivation of the scattering amplitude in terms of BS amplitude of the initial and 
final 2N states and Mandelstam electromagnetic (EM) vertex, which comprises the one-body 
and two-body parts, is proposed. Although the numerical analysis has not been performed. 
In paper |2|J] it is proposed a framework based on the BS equation approach. This is applied 
to both elastic and inelastic electron scattering. But complete calculations are performed 
within a QP framework. 

The aim of this paper is to apply the fully relativistic analysis of deuteron photodisin- 
tegration in the framework BS formalism, to segregate and estimate contribution of various 
relativistic effects in the differential cross section. 

This paper is organized as follows. In Sec. |H| we briefly discuss a connection between the 
BS formalism and QP approach, and equal time (ET) approximation. All the relativistic 
formulations of the 2N dynamics, exploiting relativistic separable interactions, are applied 
to the study of the deuteron and its inelastic observables. Basic formulae for definition 
of the Minkowski-space BS amplitude for bound and 2N continuum states are given. We 
also introduce formulae for the relativistic separable interaction kernel. Solving the BS 
equation with the separable interaction, we find the vertex function used in computation of 
the unpolarized cross section. 

In Sec. [TIT] we describe the procedure to derive deuteron photodisintegration cross section. 
The EM interaction with 2N system in the framework of the BS formalism is determined 
by the Mandelstam vertex, which generally depends on properties of the interaction kernel. 
The scheme incorporates two-body part of the Mandelstam vertex in order to guarantee the 
gauge independence. 

Sec. [IV] deals with derivation of the simplest contribution to the EM current matrix 
elements of deuteron break-up, namely the plane wave one-body approximation (PWOA). 
We discuss the transformation properties of the deuteron vertex function between the lab- 
oratory and c. m. frames. Generally one-body part of the Mandelstam vertex involves half 
off-mass-shell 7NN form factors. However, as consequence of gauge invariance the real pho- 
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ton scattering amplitude does not contain off-shell effects. In the PWOA the EM current 
matrix elements are proportional to deuteron vertex function which is taken at certain value 
of the relative energy (relative time) and 3-momentum. These are directly related to the 
photon energy-momentum transfer. Finally, we write down the expression for the differential 
photo-disintegration cross section. 

Sec. [V| is devoted to analysis of relativistic effects required by the principles of rela- 
tivity. These are relativistic kinematics and dynamics, the relative energy dependence (or 
retardation), the Lorentz contraction and spin precession. We compare results of our fully 
relativistic analysis with those of conventional non-relativistic (NR) models, the QP ap- 
proach and the ET approximation. 

In Sec. fvT] we conclude about our results. 



II. RELATIVISTIC DESCRIPTION OF TWO-NUCLEON SYSTEM 

Formulation of integral equations for amplitudes is commonly taken as the starting point 
for discussing relativistic scattering and bound state problems in the strongly interacting 
systems within quantum field theory. 

The off-shell T-matrix for the elastic scattering of two nucleons with the relative 4- 
momentum p, p' and the total momentum P satisfies to the inhomogeneous BS equation. 
In the momentum space this is the four-dimensional (4-D) integral equation with respect to 
the relative momentum k = (ko, k)Q 

T(p, p'; P) = V(p, p'; P) + ^3 fd 4 k V(p, k; P)G (k; P)T(k, p';P), (1) 

where V is an interaction kernel obtained by summing of all irreducible 2N Feynman di- 
agrams in a given field-theoretical model of the NN interaction, and Go(k; P) — the free 
two-nucleon propagator 

Go f k . P ) = [P/2 + k + m}W [P/2-fc + mp) ^ 
+ k) 2 — m 2 + it (y — k) 2 — m? + ie 

The BS amplitude for the 2N scattering state, P 2 = s > 4m 2 , is expressed in terms of 
the half off-shell T-matrix and the propagator function G as follows 

x (k; pP) = [4tt 3 z«JW(p - A;) - G (k] P)T(k, p; P)]x (0) (P; P), p-P = 0. (3) 

where p denotes on-mass-shell relative 4-momentum and x (Pj P) is the amplitude for the 
motion of free nucleons. The second term comprises rescattering contributions. 

When the two-body system has a bound state of mass M<j, the T-matrix has a pole at 

-2 
d 



P 2 = M 2 



For sake of simplicity we omit writing spinor, p-spin and polarization quantum numbers of 
amplitudes below 
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, r(p; P)r(^ ; P) 

(P ' P ; ^ ^ p2 _ M 2 » ( 4 ) 

where T is the vertex function and F is its conjugate. The vertex function of the 2N bound 
state satisfies to the homogeneous BS equation with the same interaction kernel 

T(p;P) = ^- 3 Jd 4 kV(p,k;P)^(k;P), P 2 = M 2 d , (5) 

where the bound amplitude is defined as: 

^(k;P) = Go(k;P)T(k;P). (6) 

The Eqns. ([I]) and (|5|) are manifestly Lorentz covariant and preserve the two-body elastic 
unitarity. Moreover both equations do not discriminate between the positive and negative 
energy states, yielding transformation properties of the calculated amplitudes consistent 
with the charge conjugation and time invariance. 

Although the BS calculations are feasible, a rigorous treatment of Eqs. ([I]) and (|5]) is 
rather complicated due to the appearance of a relative energy in loop integrals and the 
presence of strong singularities in the interaction kernel. The great theoretical effort have 
been applied to obtain 3-D bound-state equations from the 4-D one. A simple method to 
obtain an approximate 3-D equation has been developed by the QP approach, in which the 
BS equation is reduced to a 3-D equation by making use of a new two-nucleon propagator 
with the internal relative energy variable restricted to a fixed value. 

But there are some shortcomings in the relativistic equations obtained from the BS equa- 
tion via the 3-D reduction. First of all one meets conceptual difficulties with the consistent 
treatment of both the 2N system and its EM interactions ||25|| . Secondly by putting particles 



on mass-shell or using the positive energy projection operators, one results in a equation 
which violate the charge conjugation and CPT symmetries. Importance of the constraint 
put by discrete symmetries is discussed in Ref. |26[ . 



An alternative choice, which have been made to study the EM interactions, is the 'instant' 
or equal time approximation f24|,|2"T|j . The instant-ET approximation to EM current operator 



is a 3-D reduction consistent with charge conjugation and unitarity. This approach has 
been applied to study of the elastic electron scattering case and the deuteron breakup in 
electrodisintegration process. 

In this paper we extend the ET choice in a systematic way to describe the process of 
deuteron photodisintegration. In the PWOA the ET approximation can be obtained simply 
by replacing the initial deuteron vertex function by the BbS vertex function. The relative 
energy variable in the one-particle propagator is prescribed by the condition that the final 



state describes on-mass-shell particles [24 



The covariant BbS prescription, P ■ k = 0, puts the relative energy equal to 

k ° = 2Po^ E i p+]l ~ E \?-^' 

where = \J m 2 + k 2 and k denotes the restricted 4-vector k. The prescription leads to 
a relativistic equation of motion for the moving 2N system. By means of the boost it is 
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transformed to the Schrodinger-type equation in the rest frame[|. 

The BbS prescription is simple only in the 2N rest frame, P( ) = (y/s, 0), where the BbS 
propagator GBbs(^; P) is given by 

G Bb s(k P(o)) = -2ni5(k )G{k; P (0) )A (1) (ki)A (2) (k 2 ), (7) 

where A^(kj) is the positive energy projection operators and Q{k] P(o)) = E k (s-TE' 2 +ie) • ^ e 
delta- function in Eq. (|7|) sets the relative energy equal to zero putting both nucleons equally 
off-mass-shell. It is an attractive feature in case of the deuteron because it treats both 
nucleons in a symmetrical way and, as consequence, it is consistent with the Pauli principle. 
The QP wave function is defined in terms of the bound state vertex function T(k;P) 



defined in terms of the positive energy Dirac spinors [ [LI] : 



V m 



where Q(k; P) = j^ym 2 — k 2 . The equation satisfied by the vertex function follows from 
the BbS reduction of the BS equation: 

r(p; P) = 2^ fd 3 kV(p, k- P)Q{k- P)f(k- P), (9) 

where V(p, k; P) is a Lorentz invariant quasipotential with all relative 4-momenta are re- 
stricted by the BbS condition. The QP wave function 0qp(k; P(o)) at the rest frame of the 

deuteron, P( ) = (Md,0), is expressed in terms of Lorentz invariant \fj^4>Q~p(k] P)- By 
means of the boost relativistic equation of motion (^) for the moving deuteron transforms 
into equation in the rest frame. This invariance property yields the NR equation for the 
wave function: 



Ml - 4g 
4m 



P(o)5 P(o)) = ^ Jd 3 kV{p {0) , k( )5 P(o))^Qp(k(o)5 ^(o)), (10) 



where 3-momenta p and k in the moving frame are mapped by the boost transformation 
to p( ) and k( ) at the rest frame, respectively, and V is the quasipotential modified by 
'minimal relativity'. 

Further the discussion concerns partial decomposition of the BS vertex function of the 
deuteron T(k;P). For definiteness at the rest frame one has (we highlight dependence on 
the spin projection): 

Im(p; P(o } ) = E 9a(po, IpI; v^)r^(- P ), (M = ±1, o), = M d . (11) 



2 A detailed and systematic exposition of the covariant QP formalism for description of the elec- 
tromagnetic (EM) properties and reactions involving the deuteron is given in Ref. ]28|. 
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Here summation index a is determined by the following quantum numbers: S = 0, 1 - 
spin, L = 0, 1,2 — angular momentum, J = 1 - total angular momentum and p-spin - 
the projection of the total energy spin of the nucleon and anti-nucleon states; p) is 

the spin-angular functions and g a is the partial amplitudes. Eight partial states contribute 
to Eq. (|i~l"D. Apart from two channels with the positive energy intermediate states, viz. 
3 Sf- 3 Df, six 'extra' states, which account for anti-nucleon degrees of freedom, come into 
play. In the spectroscopic notations 25+1 Lj these are 

3 pe 3 po 1 pe 1 po 3 c~ 3 pi— 
r li r \ i r l i r l i °l i ^1 ) 

where indexes e and o stands for the even and odd parity relative to the p-spin functions. 
The partial amplitudes 1 Pf, 3 P° are even and 1 P°, 3 Pf are odd functions in the relative 
energy variable. 

Influence of admixtures of P-states and their contribution in interference with the positive 
energy states to observables in deuteron breakup and elastic proton-deuteron backward 



scattering calculated within the BS formalism is considered in Refs. |16| , p9| , |30| . It is shown 
that NN pair EM current term in the NR approach can be constructed from the P-state in 
the deuteron BS amplitude. The separate effect of the pair current term, with respect to the 
one-body, is expected to give constructive interference followed by a destructive interference 
with increasing photon energy in the range 0.1-0.7 GeV 



In the present paper we focus on the positive energy states only. So we have two channels, 
3 S^ — 3 Di , and the corresponding vertex functions can br written in the matrix form |2( 



87rr 3 / 1++ (p;P (0 )) =A/; 2 (m + pi)(l+7o)e A (m-p 2 )£ (po, |p|;s), (12) 



VmrT° t+ { V] P (0) ) = -A/p 2 (m + pi) 

3 ^ (p • e A 



x (l + 7o) \e\ + -{pi -pi) — 2 1 ( m -P2)92(Po, \p\;s) 



where s = Mj, p 1 = (P p ,p) p 2 = (E p , — p) are on-mass-shell 4-momenta, Af 1 = 

y / 2P p (m + E p ) is the normalization factor and e\ = (0, e\) is 4-polarization vector of the 
deuteron: 

+ 1 OA* pv 

H^T = -9" v + ^Pr, e A -P = 0. (13) 

A=-l 1V1 d 

We are interested in calculating matrix elements of the EM current operator between a 
state containing free nucleons and 2N bound state. To this end we need to take into account 
change in the state amplitudes, when the BbS prescription is applied to the S-matrix element 
for this process. We can utilize the normalization condition for the bound state amplitude. 
Let us suppose that the kernel V is independent of the total 4-momentum P. Then we 
obtain 



d ' k r ( p;P) dG ^ 



r(p;P). (14) 



P 2 =Mj 



2%H ^ } ' dP 2 

In terms of the radial partial vertex functions gi this conditions takes the form 
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fl L (fc ,|k|; g ) 2 (ff k -^) 

dfco / d|k| k yy^y- ; — ^ = P L , P + -?2 = 1- (15) 



2^M,y u y 1 1 ((Mi-^ + ^-fco 2 ) 2 

We find respectively for the QP vertex function g^. 



tt 2 M^ 1 1 B k (Mj-4Bj) 2 







One can deduce out of these Eqns. that the BS and QP vertex functions are formally related 
to each other as follows 



g L (k , |k|; s) 



2 



"V^^f^' (17) 



In Eq. flT7|) the QP function (|10D is related to the Schrodinger wave function by minimal 
relativity: 

0Qp(k;P (O )) = 

The presence of strong singularities in the interaction kernel can be avoided by special 
prescriptions referred to imply analytical properties. The ladder approximation with the 
interaction kernel of the BS equation V to be a sum of one boson exchange diagrams is widely 
used in solving BS equation for the NN scattering. In this case a solution of Eqn. (|5|) (for 




example, we refer to papers |Tj|- f20fl), can be found in Euclidean space after a Wick rotation. 
However that presents a substantial obstacle when one calculates observable in terms of the 
BS amplitude. Actually an analytical continuation of the solution in the complex k -plane 
back to the real axis is a procedure which carries ambiguities and extremely laborious. 
The physical solution can be obtained via the method based on the Perturbation Theoretic 
Integral Representation of Nakanishi f33fl . The BS equation for bound states is solved in 
terms of a generalized spectral representation directly in Minkowski space [33]. But the 



approach is developed for bound states in scalar theories. 

An alternative way to solve the BS equation is to use a non-local separable interaction 



kernel p3-S^|: 



N 

V(p,p'; P) = £ X ab v a {p; P)v b {p'; P), (18) 

a,6=l 

where X a b is a symmetrical matrix. In this case the bound state vertex function and the 
scattering T-matrix is obtained in Minkowski space. They exhibit analytical properties 
determined by form factors v a (p; P). The covariant form factors includes a dependence on 
P 2 , p 2 and p ■ P. Leaving the dependence on p 2 only means a simple procedure to construct 
relativistic separable interactions to be used in the BS equation. That is done in study of 
2N states in the BS formalism with the separable form of interaction, see Refs. p6| , |38| , |39| . 
These form factors are not expected to be genuine separable approximations to a realistic 
NN interactions. 
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In this paper we use the separable kernel of rank three (N = 3 in Eqn. flT8|)) for com- 
putation of the deuteron photodisintegration cross section. This interaction kernel is a 
relativistically covariant generalization of the NR Graz-II potential for the description of 
the phase shifts of the NN scattering in the coupled 3 Si- 3 /}i waves (details can be found 
in [0). The analytical properties of the radial vertex function is determined by poles in 
the relative energy: 

9o (k , |k|; s) = Ms) { l~_^ i)2 + -/%,)» ' 
g 2 (k , |k|; s) = C(s) — _ 2 j _ 2 , s = M 2 , (19) 

where A; 2 = fcg — k 2 , the coefficients A, B C are determined by the homogeneous set of 
three algebraic equations, the parameters (3^ and j a are chosen to reproduce 3 Si- 3 Di NN 
scattering phase shifts up to a laboratory energy of i?Lab — 500 MeV, the low-energy NN 
scattering parameters and the static deuteron properties (the binding energy, the quadrupole 
and magnetic moments). We adopted the parameters of the interaction kernel corresponding 
to the value of 3 _D^-state probability P 2 = 4 % and 6 %. The actual parameters (3^ and 7 a are 
chosen to have such values that the resulting T-matrix Tll'(po, P,Pq, p'; s) satisfies the exact 
two-body unitarity relation at least up to a nucleon kinetic energies i^Lab = ( m + 
with (3\i = 231 MeV, which is about 500 MeV in the laboratory system. 

This is repeated using the QP approximation to the BS equation. The QP vertex function 
is the solution of the homogeneous BS equation with the BbS version of the Green func- 
tion (|7|) and the interaction kernel equivalent to the NR Graz-II potential. The values of the 
parameters \ a b are changed relative to the BS parameters to reproduce resulting deuteron 
properties, low-energy scattering parameters and the phase shifts. The resulting radial part 
of the deuteron vertex function gi depends on the relative 3-momentum as follows 

l + 7l k 2 ~ V 2 



(k 2 + /3 1 2 1 ) 2 v 7 (k 2 + /3 2 2 ) 2 ' 

^ Wi,)a6w ra^ ' s = Ml (20) 

The particular feature of the BbS reduction is that all propagators are reduced to the 
static form. Speaking in language of the meson-exchange model, the BbS reduction com- 
pletely ignores the retardation or the relative energy dependence in the BS amplitude which 
arises from non-instantaneous effects in the NN interaction. 



III. DEUTERON PHOTODISINTEGRATION CROSS SECTION 

Let us consider disintegration of a deuteron with total 4-momentum K by a photon 
with 4-momentum q^, q 2 = 0, into a free neutron-proton (np) pair, characterized with the 
total and relative 4-momenta P and p, respectively. In the rest frame of the np pair, i. e. 
P/q) = (-y/s, 0), p M = (0,p), where y/s — is the total energy of the pair, the differential 
absorption cross section of a photon with energy ui can be written as 
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with a 



7(4*) 



the transition matrix element 



is the fine structure constant, M.'ji — the invariant amplitude, which is 

(/ I I i) °f t ne EM current operator between 
the deuteron bound state and the 2N continuum; is a photon polarization 4-vector with 
A = ±1. Momentum conservation at the photon-deuteron vertex gives K + q = P. 

Since polarizations of the particles involved in the process will not be considered here, 
averaging and summing over the photon and nuclear polarizations in the initial and final 
states, respectively, are assumed. We may choose such a coordinate system where the photon 
3-momentum is along a Z axis: q^ = (u,0,0,u). In the laboratory system, being the rest 
frame of the deuteron, the deuteron 4-momentum K/q) = (M d , 0) and the photon energy is 
denoted as £7 7 . 

In experiments on two-body photodisintegration of the deuteron the differential cross 
section ( pT|) is viewed as a function of the lab photon energy E y and angle P between 
incoming-photon and outgoing-proton 3-momenta in the c. m. system of the np pair. One 
can obtain the following kinematic relations between the laboratory photon energy and 
variables in the c. m. frame: 



iPl 



— — rrr 



s = Mi + 2E^M d 



M d 



7" 



(22) 



Following the Ref. |23| the invariant amplitude .M^ can be written in terms of the BS 
amplitude of the initial (|B|) and final states as follows: 



47T 3 



<rkd*lxsm.(hpP)W(l, k; P, K)^ M (k; K), 



(23) 



where S = 0, 1 is the total spin of the np pair and m s is its projection on to the Z axis, M 
is a projection the total angular momentum of the deuteron; A M denotes the Mandelstam 
vertex which determines the EM interaction with 2N system in the framework of the BS 
formalism. 

Let us make some remarks on current conservation. The Mandelstam vertex consists of 
one- and two-body parts, A^(p, k; P, K) = AW(p, k; P, K) + A^(p, k; P, K). The second part 
of the Mandelstam vertex determines two-body contributions of the conserved EM current. 
The specific form of A|^(p, k; P, K) depends on a given model for the interaction kernel in 
the BS equation and it cannot be associated with the pair and meson exchange currents 
in the NR approach. The gauge independence of EM current transition matrix element, 
q ■ Mfi = 0, will be fulfilled if the Mandelstam current meets the following relations: 



iq 



A.W(p,k',P,K) 



7T p (l)5 [p-k-- 



P 



-1 



+ n p (2)5[p-k + 



-l 
k) - 



s<2) (j 



P 



-i 



s(2) ( Y~ k ^ 1 



(24) 



(25) 



and 
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iq- A^(p,k;P,K) = 

E k(Ovfp + (-i 

2=1,2 L V 



,i9 



k;K)\-V[p,k- 



-i) , f;#>p(Q 



(26) 



where S^ l \p) is the fermion propagator and 7r p (Z) = 1/2[1 + r z (/)] is the projector on to 
proton state. Moreover the BS amplitudes for the initial and final states have to satisfy to 
the BS equations with the same interaction kernel. 

In Ref. [[36|] it was shown that the gauge independence condition for the elastic electron- 
deuteron scattering amplitude is fulfilled in the impulse approximation for many-rank sep- 
arable BS kernels, implying that q ■ Mjj = 0. In case of the deuteron breakup in the final 
2N state isospin 1=1 states are present, yielding a nonvanishing isovector contribution, 
q ■ M-fj 7^ 0. Consequently the two-body Mandelstam current operator should be added to 
guarantee the gauge independence. 



IV. PLANE WAVE ONE-BODY APPROXIMATION 

The amplitude (|23|) contains the FSI contributions of the final np pair. This means 
that the half off mass shell NN scattering T-matrix is needed, v. s. Eqn. (|3|). In this paper 
rescattering contributions and the pair processes in the framework of the BS formalism will 
be not taken into account. The neglect of FSI is a shortcoming of our present work. At 
present this is bound up with computational difficulties. In a forthcoming paper will be 
considered the FSI interactions from J = 0, 1 S'o, and J = 1, 3 S\- 3 Di, channels and work 
on the pair processes is in progress. Full and detailed analysis of the FSI cannot clearly be 
avoided without reconsidering an entirely different interaction kernel. 

In this investigation we confine ourselves to the simplest contribution to the EM current, 
where the photon couples to one of the two nucleons in the deuteron and the FSI between 
the outgoing nucleons are dropped — the plane-wave one-body approximation. 



A. The BS amplitude for the continuous spectrum 



According to the approximations we made the BS amplitude of the final state given in 
Eqn. (0) is the antisymmetric combination of two free Dirac positive energy spinors Q: 



Xsm s (&o,k; v^p) = 4tt 3 5(/c ) 



(27) 



XSr, 



where Xsm a (p) 



,(p)fa> + Vi)S {3) (k - p) + (-l)^ XSms (- p )(r) - Vl )S^(k + p) 



C?r s i\ u x v (p)u\ n (— p); Vo and m stands for isopin singlet and 



1 2"P2 



triplet functions respectively. Since the outgoing nucleons are on mass shell, we have con- 
straint P ■ p = which keeps the relative energy p to be equal to zero in the rest frame of 
the np pair. 



3 We use the covariant normalization of the Dirac spinors, u + u = 2E 
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B. The BS amplitude for the bound state 



The BS equation for the deuteron is solved in its rest frame. Since the vertex function 
in Eqn. (E3|) is referred to a moving frame, it has to be boosted to its rest frame. The 



Lorentz transformation between the laboratory and c. m. frames is given by 



ftx pv 



where K, 



(o) 



and Pfa 



4-momenta of the deuteron and np-pair in their rest 



frame, respectively. As only boost along the Z axis is needed, an explicit expression for the 
matrix C?^ is given by 



/ VTTv -^77 ^ 

10 

1 

V -y/rj V^+VJ 



(28) 



where the Lorentz-factor 7 = a/1 + 77 and v 7 = ^/fj, with v being the velocity of c. m. frame 
in the laboratory frame are defined in terms of the dimensionless boost parameter 77: 



v 7 ^ 



1 + 77 



E y + M d 



Under the boost the vertex function is transformed according to the general rule of the 
transformation of spinor amplitudes 



T M (k;K)=A(C)T M (C- 1 k;K l 



(0), 



(29) 



with A(£) = A^(C)A^ 2 \C), where A®(C) is the boost operator in the spinor space of Zth 
nucleon corresponding to C: 



A (/) (£) 



1 + y/T+rj" 



7o73^ 



(0 



1 + VTTv, 



(30) 



At 77 — > the matrix fl28|) and operator (|30|) turn to unity, i. e. C — > I and A —> I. This 
corresponds to the static limit for the BS amplitude and takes place on threshold of deuteron 
photodisintegration, of ^ < 1 (v. s. Eqn. (g|)). In the Eq. © the 7073-term affects 

the spin degrees of freedom of the BS amplitude. 



C. The electromagnetic vertex 

One-body part of the Mandelstam vertex has the form 

a w ( P , k;P,K) = <«yw (p - k - § ) r« (| + P , *L + k) is® (I - k)' 1 

+ «w ( P - fc + 1) (f - P , I - fc) (| + fc ) " , (3i) 

where r^'(p, k) is off-mass-shell 7NN vertex for the Zth nucleon. This is common problem 
of direct application of the BS equation to the NN interaction. The consequence of gauge 
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constraints for off-shellness in the EM vertices have been recently considered in Ref. f40 
(see also references therein). In our case we deal with half off- mass shell EM vertex. This 
type of vertex occurs in the (e, e'N) reaction, e. g. nucleon knockout and inclusive electron 
scattering, when the initial nucleon is taken to be bound (off-mass-shell) and the knocked- 
out nucleon is assumed to be in physical state. It is shown in Ref. [[|0| that the off-shell 
behavior of the EM vertex of the nucleon in a real Compton scattering on a free nucleon 
does not play a role as consequence of gauge invariance. 

Thus we deal with well-known on-mass-shell version of the EM vertex: 

r<?(?) = 7, {Fl s \q) + t 3 ( V(%)) + + r®Fl%)) , (32) 

where Fi^ (q) is the isoscalar and isovector Pauli-Dirac form factors of the Ith nucleon, nor- 
malized as: Fi s) (0) = §, F 2 (s) (0) = ^±^, F^O) = §, F 2 M (0) = with the anomalous 
part of the proton and neutron magnetic momenta respectively. 



Employing the written transformation laws, substituting Eqns. (|27|) , fl29|) , fl31|) to 
Eqn. ( P^) and integrating over intermediate 4-momentum, we arrive at an expression in 
terms of the c. m. vertex function and propagators 

M£ = £ xil(0,p;v^p)r«(g 2 = 0)A(/:)^(fc;K )r M (A; 0i ,k / ;K ), (33) 
1=1,2 

where k\ = C~ x {p + (— 1)'|). All possible contributions to the transition matrix element in 
the PWOA are depicted in Fig. [TJ. Since both the initial and final 2N states are antisymmetric 
the four diagrams are identical. 

In the PWOA it is seen that the matrix element is proportional directly to the deuteron 
vertex function taken at specific value of energy-momentum. As the relative 4-momentum 
ki is restricted by energy-momentum conservation in a photon-nucleon vertex, the relative 
energy k$i and 3-momentum k; variables depend on the photon energy in the lab. frame 
and the dimensionless boost parameter rj: 

koi = Vv\p\\\ + (-l) 1 ^, 

k± z = p±, u = M d ^], (34) 

m<1 



k]|/= v/l + r/P|| + (-l) Y 

where indices 1 1 and _L denotes the longitudinal and transverse components of vector k with 
respect to direction of the incoming photon 3-momentum q, here |q| = E 7 . The situation 
is illustrated on Fig. 2. For photon energies E 1 < 0.2 GeV one is probing the energy- 
momentum distribution of the bound nucleons in the deuteron where the high-momentum 
'tails' of the nucleonic states in deuteron are especially relevant. For given c. m. angles O p 
it is found that both the relative energy, which accounts for the retardation in the vertex 
function of the deuteron, and the modulus of the 3-momentum of the proton (neutron) 
rise strongly with E 1 . They are smaller at forward scattering angles, and for other angles 
covering wide bands from 100 MeV to 2 GeV and from 500 MeV to 2.5 GeV. 
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Introducing the deuteron state components, the 2N continuum amplitude in the matrix 
representation (for details we refer to the paper |20|), the transition matrix element can be 
evaluated calculating traces of 7-matrix expressions: 



M£ = - Sp( x Sms (p)r PiM A(£)^(-^ + fc l5 K )T M (k i; K )A(£- 



Sp ( Xs ms (p)A(£)r A/ (fc 2 ; K )S(^- - k 2 ; K^AiC- 1 )^ ) , (35) 



with 



M 2 

XimAv) = ^(^-P2)C s ( 1 + 7o)(^ + Pi), (36) 

K 

Xoo(p) = — ~^( m ~ ^2)75(1 + lo)(m + p 1 ), 

where £ ms is a polarization 4-vector with the following completeness and orthogonality re- 
lations 

E c s c; = -r + ^ e-p = o, (37) 

m s =— 1 

and normalization constants M v and vectors pi^ are defined in Eqn. flT2l). 

Since the EM nucleon form factors can be taken by their on-shell form, we have for the 
charge-current operator 

r PlM = 7^ + 2^<V9 » F «,m = 2^<V9 • ( 38 ) 
The fermion propagator in Eqn. fl55]) 

= p^T r ( 39 ) 
k z — m z + xe 

is connected with the propagator S by S = —CS T C, where C = ry 2 7°. 

Using the expression (P5Q, we find that the differential photo-absorption cross section 
can be written as 

da da dcr SP 



d/^^iip d^\<p cii^^t'p 

where the <7o is the part of the cross section which makes up the shape of the angular 
distributions 



da a ( 1 + VI + 2 
dQ p 47rs \ 2 y s _ n , 



E I Xo S I 2 (41) 



and the rfcxsp accounts for the effect of the boost on the spin degrees of freedom (the spin 
precession) of the nucleons 



14 



dasp _ a ( 1 + VI + ' 
dVtp A-ks \ 2 



with the square modulus of the amplitude A^fp is given by 



/ 5=0,1 



S |2 
SP I 



(42) 



| A^fp | = 2/3Re(X^X^*) 

+ [3 2 (| Xf | 2 -2 J Re(X 5 Xf)) - 2/3 3 J Re(XfX 2 5 *) + /3 4 | X 2 5 | 2 (43) 

with (3 = 1+ ^ +T? and the amplitudes Xf (i = 0, 1, 2) are expressed as 



X 5 = ^^5 ms (p)r p , A ^ 1 )( Sl ;iro)r A f(A;i;Xo)J + n, 
Xf = Sp(xsm s (p) r p>A 7on 3 S' (1) (si;Xo)r M (A; 1 ;Xo) 

- 5 , p^X5 ms (p)r p ,A5' (1) (- s i;^o)rAf(^i;Xo)7on3j + p^ n, (44) 
X 2 5 = Sp[xsm s {p)r pA 7onzS {1 \s 1 ]K )r M {ki,K )'y n 3 J + p 



n 

with n 3 = (0, 0, 0, 1) denoting a unit vector and si = K /2 — k%. 

These are the general expressions for the deuteron photodisintegration cross section (21) 
in the PWOA. Since we omit the two-body contribution to the transition matrix element, 
we does not preserve the gauge independence of the amplitude. Thus averaging over photon 
polarizations, we make use of Coulomb gauge, e° = 0, e ■ q = with the completeness 
relation of the form 

E ( £ x)i(£x)j = <% _ ^T"> hj=x,y. (45) 

A=±l Q 



All 7-matrix expressions in the matrix elements in Eq. ( 44]) and the square modulus of 
the amplitude in Eq. (|43|) are evaluated with the computer algebraic program REDUCE. 
When summing over nuclear polarizations, we make use of the relations ( |T3"D and ([37]). 



V. ANALYSIS OF RELATIVISTIC EFFECTS 

Now we are in a position to do the final computations. The results for the angular 
distributions for deuteron photodisintegration at four different photon energies E y in the 
laboratory frame are depicted in Fig. |3|. The cross section is calculated in the framework 
of the BS equation with the separable interaction kernel (with two different strength of D- 
state, P2 = 4 % and P2 — 6 %). This interaction kernel is similar to that employed in the 



calculations of the deuteron EM elastic form factors ||36|| . For simplicity, though it might be 
important, in our calculations we disregard negative-energy partial states of the deuteron 
vertex function. These will be considered in detail in the BS formalism for a meson exchange 
interaction kernel. 

In Fig. ^| it is seen that the resultant curves reproduce shapes of angular distributions [|TJ . 
Above the threshold it is an almost perfect sin 2 B p behavior, which corresponds to El 
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transition to the 3 Pq np continuum state. At higher photon energies the overall magnitude 
of the cross section rapidly falls off and the maximum is shifted from 90° to 70° and, further, 
to 60°. The whole distributions, which are dominated by magnetic transitions, become 
flatten, the ratio of the forward cross section to the maximum decreases. The role of the in- 
state becomes more pronounced. But notwithstanding these similarities, the theory appears 
systematically less than the experimental distributions. At E^=20 MeV it is less by factor 
1.4, at _E 7 =100 MeV by factor 3, in the A-resonance region it is expectedly lower by almost 
factor 10, and at E y =500 MeV by factor 3 with respect to .D-state weight P 2 = 6 %. In the 
NR approach, incorporating meson-nucleon degrees of freedom, a good deal of discrepancy 
between the experiment and theory is diminished by contributions from meson-exchange 
currents and isobar configurations. 



A. Approximate calculations 

Despite the fact that our theoretical results does not describe the experimental data, 
we are able to make definitive statements concerning the relative importance of the various 
relativistic effects. We distinguish the following classes of the relativistic effects in our consis- 
tent treatment: relativistic kinematics and EM current operator; relativistic NN dynamics, 
which forbids instantaneous interactions, and leads to the relative energy dependence of 
the deuteron vertex function (retardation), the boost transformation affecting the internal 
variables of the BS amplitude of the deuteron (Lorentz contraction), the part of the cross 
section denoted as ^Q, and its spin degrees of freedom (spin precession), the cross section 

dn ' 

We start discussion with evaluation of the size of contributions due to retardation, 
Lorentz contraction and spin precession. Let us consider a number of approximate cal- 
culations with respect to the exact positive-energy BS calculations: 

1. First of all we perform the static approximation (BS-SA) to the BS cross sections, 
which amounts to neglecting the boost on the arguments of the deuteron vertex func- 
tion and one-particle propagator, see Eq. fl3"3]). It is achieved by putting the boost 
parameter r\ = in the deuteron vertex function TL(k i,ki)\ rj=0 = Tl(p i,pi), where 

Poi = {—) l ^2 an d Pi = P + ( — ) Z f = 1)2). Kinematically this effect is shown in 
Fig. 0. The booster in Eq. ( PU|) is approximated by A(C)\ V=Q = I. This approximation 
excludes contributions due to the Lorentz contraction and spin precession. 

2. Moreover we pay special attention to investigation of the influence of the boost on the 



nucleon relative energy in one-particle propagator. As it is shown in Ref. |35j that this 
is the most important relativistic contributions to the deuteron EM form factors. Here 
we consider the case (BS-BR), which is the same as BS-SA, but includes the boost on 
one-particle propagator S^(k i, |p z |) due to recoil, k oi = w^|p| + 



Finally, in order to find out the relativistic correction associated with the relative en- 
ergy dependence in the matrix elements, which is brought by the BS vertex function. 
We consider the zero-order approximation (BS-ZO) for the vertex function, i. e. com- 
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puting the radial parts of the vertex function in the BS-BR approximation with koi 
equal to zero. 

In Fig. f| we present the results for the angular distributions at the same lab. photon 
energies as in Fig. |3|. The solid curve corresponds to the exact calculation. We would like to 
stress the importance of the effects concerning the booster. In Fig. [| it is clearly seen the role 
of the dasp cross section (3-dot-dash line) with rise of the energy. It becomes pronounced 
at the scattering in the forward semi-sphere, where it is almost the half of size of the cross 
section. One should allow for such a contribution starting at the medium photon energies. 
When added to the BS-BR approximation (dash line), the latter becomes to be plausible 
for discussed _E 7 region. The BS-BR approximation, which includes the boost on the one- 
particle propagator, gives the shape of the angular distributions close to the exact ones. We 
can conclude that the BS-BR approximation supplemented with dcr SP contribution accounts 
for major relativistic effects in the cross section. On the other hand the SA approximation 
(dot line) ceases to be reasonable for the E 1 above 100 MeV, as it has a wrong position of 
cross section maxima. 

We show the influence of the boost transformation in the arguments of the initial-state 
vertex function. It is worthwhile to mention that the boost leaves the arguments of the radial 
part of the vertex function unchanged (the radial function <7z,(p; K) depends on the Lorentz 
invariants K 2 , p 2 and p ■ K). It only has a direct bearing on its spin-orbital part. In Fig. |5| 
it is displayed the relative deviation of the BS-BR (dash line) and BS-ZO (dot-dash line) 
approximations from the BS cross section, da^. The deviation of the BS-BR cross section is 
due to the boost effect on the orbital part of the deuteron vertex function. The respective 
contribution is quite great especially at forward and backward proton c. m. angles. The 
effects of retardation is responsible for discrepancy between the two curves in Fig. [5[ One 
can see that this is practically uniform difference reaching up to 5 %. If one ignores the 
dependence on the relative energy in the deuteron vertex function, one comes up with the 
cross section which is slightly smaller but keeps the same shape. 



B. Comparison with other approaches 

The BS formalism consistently accounts for the relativistic effects associated with mani- 
fest Lorentz covariance of scattering amplitudes. As far as matter of relativity is concerned, 
we compare the exact results of the BS framework with the following approaches. 

1. In the PWOA the ET approximation can be obtained immediately from the exact 
expressions replacing the BS deuteron bound state for "++"-channels by the QP vertex 
function, which is solution of the 3-D QP equation with the BbS propagator (0) and 



with the refitted version of the separable interaction kernel Graz-II ||36|| . Essentially 
this approximation makes use of instantaneous interactions, i. e. with the zero relative 
time, or respectively, relative energy in the deuteron vertex function and the one- 
particle propagator. 

2. A minimally relativistic approach employs the same QP vertex function of the deuteron 
corresponding to solution of the bound state equation the BbS propagator. This 
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corresponds to the static limit for the ET approximation (ET-SA). The approach 
incorporates the relativistic kinematics and covariant form of the EM current matrix 
elements. 

3. A purely non-relativistic approach makes use wave function of the deuteron, which is 
solution of the Schrdinger equation with NR Graz-II separable potential. 

In Fig. |H we compare the above mentioned approaches with the exact relativistic calcu- 
lation. One can see that NR approach (dot line) is very crude approximation for the photon 
energies greater that 100 MeV. The minimally relativistic approach (dot-dash line) improves 
the situation. By no means the NR approach is reasonable to describe deuteron photodis- 
integration at high photon energies above the pion threshold. At least one should include 
minimally relativistic corrections. The conclusion is in accordance with the discussion in 
Ref. 1], where it was shown that should the relativistic effects not included, a theory gave 
too much peaking of the differential cross section at P =O° and 180°. 

As one can see the overall sign of the relativistic contribution to the NR cross section is 
negative, but it depends on given photon energy in the range G p =60-hl20°. 

Taking into account of Lorentz deformation (dash line in Fig. |6|) produces sizeable effect. 
It is seen that the ET approximation is a good approximation at low photon energies. Almost 
equally at all proton c. m. angles, the discrepancy between the exact and ET calculations 
reaches about 10 % at £ 7 =300 MeV and 25 percent at the £ 7 =500 MeV. 



C. Expansion of the relativistic model 

It is well-known that the usual way to include relativistic effects is — - and —-expansion 

of the exact relativistic model, if one confines to the lowest order correction and — 
beyond the NR amplitudes |lj. This relativistic correction is valid at the photon energies 
Ej up to few hundred MeV and it seems to correspond to corrections of spin-orbit type to 
the NR current operator. 

We analytically performed approximation of the expression (|43"1) for the angular distri- 
butions in the limit E 1 <C m. For the c. m. frame variables in the NR limit, corresponding 
to energies E 1 < 100 MeV, we have |p| = \jm{E 1 — e^), uj = E 7 and the boost parameter 
./fj = 0. Thus the matrix element is expanded in powers of — and — with — « 0.3 keeping 

V' ^ 1 m m m 1 

the lowest order terms, recoil effects and boost effects are neglected as well. The result is 
that of the non-covariant description of deuteron in the framework of the conventional NR 
models incorporating relativistic effects in a — expansion of a relativistic model PU3l|. 
The differential cross section of a photon with energy E 7 in the NR framework is given 

by 

^-^i^l' (46) 

where the deuteron break-up matrix elements are expressed in terms of the NR deuteron 
wave function \1/ M (k), the wave function of the 2N scattering states ^ p s ms (k) in the following 
form 
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-E 



=1,2 



r/k 

* pSms (k)^(q)* M (kO ) (A = ±l), 



(27T) 



(47) 



where kj = p — (— 1)'| and the EM nucleoli form factor is given by 

F z A (q) = (-l)('+i)i± r * pA ■ Ks + T * )Kv z[<t( ° X q] 



m 



2m 



with n s = \{k p + K n ) and /«„ = |(k p — k„). 

In the PWOA one can reduce square of the amplitude (fl"7|) in Eq. ([46|) to 'textbook' 
formulae which reproduces the shape of a angular distribution but misses its absolute size: 



El 



dao = sin2 e ^ 2 + w ") + 3^1 ^ + ^ 2 ) + ^(^ + w *) 



+ 



2m 2 



2U 1 U 2 



U X W 2 + f/ 2 W: 



i (l + 3 cos(29 p )) + ^-^(5 + 3 cos(26 p ))] } 



V2 



(48) 



where indexes 1 and 2 at the S- and D-state of the deuteron, denoted as U and W, respec- 
tively, means that they are evaluated at modulus of the proton and neutron 3-momentum 
in the deuteron, k x and k 2 . These are formally related to the BS form factors as 



U 



'irm 



2E k - M d 



g (0,\k\;s = M ( 



2 ) 



W 



^7im 
2£ k - M d 



g 2 (0,\k\;s = M 2 d ). 



VI. CONCLUDING REMARKS 

The objective of the present study is evaluate the various relativistic contributions to the 
angular distributions in deuteron photodisintegration process. That can be done in a variety 
of theoretical frameworks and dynamics. In this paper we have applied the fully relativis- 
tic formalism, based on the Bethe-Salpeter equation for the 2N scattering amplitude and 
deuteron bound state. Beyond the choice of the theoretical framework, which is manifestly 
covariant at every step of the calculation, the important issue is dynamical model of the 
nucleon-nucleon interaction. Pursuing the aim to obtain clear understanding and conduct 
straightforward comparison with the non-relativistic and minimally relativistic approaches, 
we employ the effective separable interactions in construction the solvable dynamical model 
of the deuteron. 

In order to obtain ultimate results in analytical form we discard channels containing 
negative-energy states in the Bethe-Salpeter amplitude of the deuteron. We also do not 
include the two-body contributions to the electromagnetic current operator and neglect the 
final state interaction in the outgoing 2N state. The last two limitations constitutes the so- 
called plane wave one-body approximation. The negative-energy states, or P-states, in the 
deuteron vertex function are not presumably irrelevant. Their contributions are expected 
to be significant within the considered interval of the photon energies. In this respect our 
present study is primarily of a comparative character. 
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Despite to these specifications, the strongest advantage in our investigation concerns the 
fully covariant and rigorous description of the bound state and the deuteron electromagnetic 
current. The present approach accounts for wealth of the relativistic effects to the differ- 
ential cross section of deuteron photodisintegration: the role of relativity of the transition 
matrix elements between nuclear states, the influence of the retardation in the deuteron 
vertex function and one-particle propagator and changes in amplitudes due to the Lorentz 
deformation and spin precession. 

Incorporation of these relativistic effects can play a crucial role in theoretical analysis of 
deuteron photodisintegration even at intermediate laboratory photon energies for the for- 
ward and backward scattering. Here the most important contributions comes from the boost 
in the arguments of the initial state vertex function and the boost on the relative energy in 
the one-particle propagator due to recoil. As one is concerned with the covariant approaches, 
the equal time approximation is more or less reasonable approach from a pragmatic point 
of view. 

Further, we can draw the following conclusions out of the present investigation: 1) the 
Bethe-Salpeter approach allows one to take into account Lorentz invariance and relativistic 
dynamical structure of the two-nucleon system in the most general form. 2) The novel 
feature brought by the Bethe-Salpeter approach is the retardation due to the dependence 
of the bound state amplitude on the relative energy of the nucleons. In the plane wave 
one-body approximation, the scattering amplitude bears the explicit dependence on this 
variable and its magnitude is measured by the photon energy. In our opinion, this turns out 
to be by far the most important fact enabling to study recoil effects due to energy transfer 
to a nucleon by a photon. 3) The role of the boost transformation of the spin degrees of 
freedom becomes noticeable in increasing order of the photon energy at forward scattering. 

Finally, the region of high photon energies (above i? 7 =500 MeV) calls for a more complete 
investigation. In this energy region one needs to construct a realistic interaction kernel in the 
Bethe-Salpeter equation. Moreover, extending the above calculations includes contributions 
due to P-states and the two-body processes in the EM current operator matrix elements. 

ACKNOWLEDGMENTS 

We are indebted to V.V. Burov, S.G. Bondarenko, A.A. Goy, A.V. Molochkov and A.V. 
Shebeko for many useful discussions and encouragement. This work has been partially 
supported by grant 02.01.22/2000, "Universities of Russia". 



20 



REFERENCES 



[1] H. Arenhovel and M. Sanzone, Few-Body Systems Suppl. 3, (1991). 
[2] D. Babusci et al. Nucl. Phys. A 633, 683 (1998). 

[3] F. Adamain, Abstracts of Xlth Particle and Nuclei Int. conference. Uppsawa. Sweden. 
(1999) 93. 

[4] J.E. Belz, Phys. Rev. Lett. 74, 646 (1995). 

[5] H. Ejiri, Private communication. XIV Int. Seminar on High Energy Physics Problems. 
Dubna. 1998. 



[6] C. Bocha et al, |nucl-ex/980800Tl 11 Aug. 1998. 
[7] R. Crawford et al, Nucl. Phys. A 603, 303 (1998). 
[8] S. Wartenberg et al., Few-Body Systems 26, 213 (1999). 
[9] P. Wilhelm and H. Arenhovel, Phys. Lett. B 318, 410 (1993). 
[10] J.M. Laget, Nucl. Phys. A 579, 333 (1994). 

[11] W. Jaus, B. Bofmger and W.S. Woolock, Nucl. Phys. A 562, 477 (1993). 
[12] J.M. Zuilhof and J.A. Tjon, Phys. Rev. C 22, 2369 (1980). 
[13] A. A. Logunov and A.N. Tavkhelidze, Nuovo Cim. 29, 370 (1963). 
[14] R. Blankenbecler and R. Sugar, Phys. Rev. 142, 1051 (1966). 
[15] J. Adam, Jr, J.W. Van Orden and F. Gross, Nucl. Phys. A 640, 391 (1998). 
[16] L.P Kaptari, B. Kampfer, S.M. Dorkin and S.S. Semikh, Phys. Rev. C 57, 1097 (1998). 
[17] A.Yu. Umnikov, F.C. Khanna, K.Yu. Kazakov and L.P. Kaptari, Phys. Lett. B 334, 
163 (1994). 

[18] C. Ciofi degli Atti, L.P Kaptari and S. Scopetta, Eur. Phys. J. A 5, 191 (1999). 

[19] N. Honzawa and S. Ishida, Phys. Rev. C 45, 17 (1992). 

[20] L.P. Kaptari et al.., Phys. Rev. C 54, 986 (1996). 

[21] A.V. Anisovich and V.A. Sadovnikova, Eur. Phys. J. A 2, 199 (1998). 

[22] A.E.L. Dieperink and S.I. Nagorny, Phys. Lett. B 456, 9 (1999). 

[23] A.Yu. Korchin and A.V. Shebeko, Sov. J. Nucl. Phys. 54, 375 (1991). 

[24] E. Hummel and J.A. Tjon, Phys. Rev. C 49, 21 (1994). 

[25] D.R. Phillips and S.J. Wallace, Few-Body Systems 24, 175 (1998). 

[26] V. Pascalutsa and J.A. Tjon, Few-Body Systems Suppl. 10, 105 (2000). 

[27] D.R. Phillips, S.J. Wallace and N.K. Devine, |arXiv:nucl-th/9906068| . 

[28] W. Jaus and W.S. Woolock, Helv. Phys. Act. 57, 644 (1984). 

[29] S.G. Bondarenko, V.V. Burov, M. Beyer and S.M. Dorkin, Phys. Rev. C 58, 3143 (1998). 
[30] L.P. Kaptari, B. Kampfer, S.M. Dorkin and S.S. Semikh, Few-Body Systems 27, 189 
(2000). 

[31] M. Anastasio and M. Chemtob, Nucl. Phys. A 364, 219 (1981). 
[32] G. Rupp, Nucl. Phys. A 508, 131 (1990). 

[33] N. Nakanishi, Graph Theory and Feynman integral, Gordon and Breach, New York, 
1971. 

[34] K. Kusaka, K.M. Simpson and A.G. Williams, Phys. Rev. D 56, 5071 (1998). 
[35] G. Rupp and J.A. Tjon, Phys. Rev. C 41, 472 (1990). 
[36] G. Rupp and J.A. Tjon, Phys. Rev. C 45, 2133 (1992). 

[37] K. Schwarz, J. Haidenbauer and J. Frohlich J., Phys. Rev. C 33, 456 (1986). 
[38] G. Rupp and J.A. Tjon, Phys. Rev. C 37, 1729 (1988). 



21 



[39] S.G. Bondarenko, V.V. Burov and S.M. Dorkin, Proc. XII Int. Seminar on High Energy 

Physics Problems/Ed Baldin A.M., Burov V.V. Dubna: JINR (1994) 227. 
[40] S.I. Nagorny and A.E.L. Dieperink, Eur. Phys. J. A 5, 417 (1999). 



22 



FIGURE CAPTION 



FIG. 1. Diagrams corresponding to the plane wave one-body approximation to the matrix 
element of deuteron photodisintegration. Outgoing particles are on their mass-shell. 

FIG. 2. Modulus of the energy \ko\ and 3-momentum |k| of a nucleon in the deuteron 
in its rest frame versus the photon energy E 1 (solid lines). Dash lines correspond to solid 
ones while neglecting the boost parameter, i. e. ^Jfj = (static approximation). A set 
of fixed proton scattering angles O p = 0°, 90 and 180° in the c. m. frame correspond to 
curves labeled as 1, 2 and 3, respectively. The same curves are related to the energy and 
3-momentum of a spectator nucleon at cm. frame angles O p = 180°, 90 and 0°. 

FIG. 3. The differential cross section in the plane wave one-body approximation at 
different photon energies Ej. Curves corresponds to different probabilities of 3 Di partial 
state. Solid line, P2 = 4 %, dash line, P2 = 6 %. 

FIG. 4. The differential cross section in the plane wave one-body approximation at 
different photon energies E 1 . Curves: solid line — the exact positive-energy BS calculation, 
dotted line — the static approximation (exclusion of the Lorentz contraction), dash line 
- the static approximation with taking into account the boost on one-particle propagator 
due to recoil (the Lorentz contraction), 3-dot-dash line — the contribution due to the spin 
precession, Eqn. (142]). Probability of 3 Df partial state is P2 = 4 %. 

FIG. 5. The relative deviation of the deuteron photodisintegration cross section for 
the photon energies E 1 = 100, 300 and 500 MeV for the approximations with respect to 
the BS result: BS-ZO (dash line) and BS-BR (dot-dash line). On the Y-axis it is plotted 
ff °~^ o BS -" x 100 %, where a =BR, ZO. Probability of 3 Df partial state is P 2 = 4 %. 

FIG. 6. The relative difference of the deuteron photodisintegration cross section for 
the photon energies E 1 = 100, 300 and 500 MeV for the following cases: the equal time 
approximation (dash line), the minimally relativistic approach (dot-dash line) and non- 
relativistic approaches (dotted line) with respect to the exact positive-energy BS result. On 
the Y-axis it is plotted x 100 %, where a =ET, ET-SA and NR. Probability of 3 Df 
partial state is P2 = 4 %. 
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FIG. 3. Kazakov K.Yu. 
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FIG. 4. Kazakov K.Yu. 
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